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1 Classical orthogonal polynomials 



One aspect in the theory of orthogonal polynomials is their study as special functions. 
Most important orthogonal polynomials can be written as terminating hypergeometric series 
and during the twentieth century people have been working on a classification of all such 
hypergeometric orthogonal polynomial and their characterizations. 

The very classical orthogonal polynomials are those named after Jacobi, Laguerre, and 
Hermite. In this paper we will always be considering monic polynomials, but in the literature 
one often uses a different normalization. Jacobi polynomials are (monic) polynomials of 
degree n which are orthogonal to all lower degree polynomials with respect to the weight 
^ ■ function (1— on [—1, 1], where a, (3 > —1. The change of variables x ^— 2x— 1 gives 
Jacobi polynomials on [0, 1] for the weight function w{x) = x^{l — x)", and we will denote 
these (monic) polynomials by Pn"'^\x). They are defined by the orthogonality conditions 



[ Pj^'''f^\x)xf^{l-x)''x''dx = 0, A; = 0,1,... ,n- 1. (1.1) 
Jo 



The monic Laguerre polynomials Li^\x) (with a > —1) are orthogonal on [0, oo) to all 
g ■ polynomials of degree less than n with respect to the weight w{x) = x"e~^ and hence satisfy 
the orthogonality conditions 

: [ L(,")(x)x"e-^x^da; = 0, k = 0,1, . . . ,n - 1. (1.2) 

Finally, the (monic) Hermite polynomials Hn{x) are orthogonal to all lower degree polyno- 
mials with respect to the weight function w{x) = on (— oo, oo), so that 

Hn{x)e-''\^ dx = 0. A; = 0,1,... ,n- 1. (1.3) 

) 

These three families of orthogonal polynomials can be characterized in a number of ways: 

• Their weight functions w satisfy a first order differential equation with polynomial 
coefficients 

a{x)w\x) = p{x)w{x), (1.4) 



*The first author is a Research Director of the Belgian Fund for Scientific Research (FWO). Research 
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with a of degree at most two and p of degree one. This equation is known as Pear- 
son's equation and also appears in probabihty theory, where the corresponding weights 
(densities) are known as the beta density (Jacobi), the gamma density (Laguerre), and 
the normal density (Hermite). Note however that for probability density functions 
one needs to normalize these weights appropriately. For the Jacobi weight we have 
(t(,t) = ,x(l — x), for the Laguerre weight we have (j{x) = x, and for the Hermite weight 
wc see that a{x) = 1, so that each family corresponds to a different degree of the 
polynomial a. 

The derivatives of the very classical polynomials are again orthogonal polynomials of 
the same family but with different parameters (Sonin 1887, W. Hahn 1949). Indeed, 
integration by parts of the orthogonality relations and the use of Pearson's equation 
show that 

^Hn{x) = nHn_i{x). 
ax 

The differential operator D = d/dx therefore acts as a lowering operator that lowers 
the degree of the polynomial. 

Pearson's equation also gives rise to a raising operator that raises the degree of the 
polynomials. Indeed, integration by parts shows that 

^[x^{l-xrPi''^^\x)\ = -^a + P + n)xP-\l-xyp!;,l\^^^-^\x), (1.5) 
£ [x"e-^4")(x)] = -x--\-^L^:;^\xl (1.6) 



d 

dx 



e-- H^{x) = -2e-- H^^^{x). (1.7) 



The raising operator is therefore of the form a{x)/w{x)Dw{x). Using this raising 
operation repeatedly gives the Rodrigues formula for these orthogonal polynomials: 

[x^+"(l - x)"+"] = (-l)'^(a + /? + n+l)„x^(l-x)"Pi"''^)(x), (1.8) 



dx 

^[a;"+"e--] = (-l)"a;«e--L(")(a;), (1.9) 

ile--' = (-l)"2-e--'//„(x). (1.10) 
dx^ 

The Rodrigues formula is therefore of the form 

— [a''{x)w{x)] = CnW{x)Pn{x), 

where C„ is a normalization constant (Hildebrandt 1931). 
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• Combining the lowering and the raising operator gives a hnear second order differential 
equation for these orthogonal polynomials, of the form 

a{x)y"{x) + T{x)y'{x) = \ny{x), (1.11) 

where a is a polynomial of degree at most 2 and r a polynomial of degree at most 
1, both independent of the degree n, and An is a constant depending on n (Bochner 
1929). 

The Laguerre polynomials and the Hermite polynomials are limiting cases of the Jacobi 
polynomials. Indeed, one has 

lim a"Pi"''^)(x/a) = lI^\x), (1.12) 

and 

lim 2"a"/2pK«)(^±^) = H^i^x). (1.13) 
a^oo 2ya; 

The Hermite polynomials are also a limit case of the Laguerre polynomials: 

lim {2a)-''/^Ll^\V2^x + a) = Hn{x). (1.14) 

a— >oo 

In this respect the Jacobi, Laguerre and Hermite polynomials are in a hierarchy, with Jacobi 
leading to Laguerre and Laguerre leading to Hermite, and with a shortcut for Jacobi leading 
to Hermite. This is just a very small piece in a large table known as Askey's table which also 
contains classical orthogonal polynomials of a discrete variable (Hahn, Meixner, Kravchuk, 
and Charlier) for which the differential operator D needs to be replaced by difference oper- 
ators A and V on a linear lattice (a lattice with constant mesh, see [^). Finally, allowing 
a quadratic lattice also gives Meixner-Pollaczek, dual Hahn, continuous Hahn, continuous 
dual Hahn, Racah, and Wilson polynomials, which are all in the Askey table. These poly- 
nomials have a number of g-extensions involving the g-difference operator and leading to 
the g-extension of the Askey table. In Andrews and Askey suggest to define the classical 
orthogonal polynomials as those polynomials that are a limiting case of the 4(y93-polynomials 

' aqMq.cq ' ^' ^ 

with A(x) = + q^^^cd and hdq = (these are the g- Racah polynomials) or the 4993- 
polynomials 



a"lV„(x; a, 6, c, (i|g) / q^",q'^^^abcd,ae'-^,ae 



{ab;q)n{ac;q)n{ad;q)n \ ab,ac,ad 

with X = cos 6 (these are the Askey- Wilson polynomials). All these classical orthogonal 
polynomials then have the following properties: 

• they have a Rodrigues formula, 

• an appropriate divided difference operator acting on them gives a set of orthogonal 
polynomials. 
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they satisfy a second order difference equation in x which is of Sturm-Liouville type. 



The classical orthogonal polynomials in this wide sense have been the subject of intensive 
research during the twentieth century. We recommend the report by Koekoek and Swarttouw 
P^ ], the book by Andrews, Askey and Roy 0, and the books by Nikiforov, Uvarov [^, and 
Nikiforov, Suslov, and Uvarov |^ for more material. Szego's book |Q is still a very good 
source for the very classical orthogonal polynomials of Jacobi, Laguerre, and Hermite. For 
characterization results one should consult a survey by Al-Salam 



2 Multiple orthogonal polynomials 

Recently, there has been a renewed interest in an extension of the notion of orthogonal 
polynomials known as multiple orthogonal polynomials. This notion comes from simulta- 
neous rational approximation, in particular from Hermite-Pade approximation of a system 
of r functions, and hence has its roots in the nineteenth century. However, only recently 
examples of multiple orthogonal polynomials appeared in the (mostly Eastern European) 
literature. In this paper we will introduce multiple orthogonal polynomials using the or- 
thogonality relations and we will only use weight functions. The extension to measures is 
straightforward. 

Suppose we are given r weight functions wi, W2, • • • , on the real line and that the 
support of each Wi is a subset of an interval Aj. We will often be using a multi- index 
n = (ni, 77,2, . . . , rir) G N*" and its length |n| = rii + n2 + ■ ■ ■ + rir- 

• The r- vector of type I multiple orthogonal polynomials (A^^i, . . . , A^,r) is such that 
each j is a polynomial of degree rij — 1 and the following orthogonality conditions 
hold: 

I x^^A,n,j{x)wj{x)dx = Q, A; = 0, 1, 2, . . . , |r?| - 2. (2.15) 

Each Afi^i has Ui coefficients so that the type I vector is completely determined if we 
can find all the \n\ unknown coefficients. The orthogonality relations ( [^.15| ) give \n\ — 1 
linear and homogeneous relations for these \n\ coefficients. If the matrix of coefficients 
has full rank, then we can determine the type I vector uniquely up to a multiplicative 
factor. 

• The type II multiple orthogonal polynomial is the polynomial of degree \n\ that 
satisfies the following orthogonality conditions 

0, A; = 0, 1,... ,ni - 1, (2.16) 
0, A; = 0,1,... ,n2-l, (2.17) 

0, A; = 0,1,... ,n^-l. (2.18) 

This gives \n\ linear and homogeneous equations for the \n\ + 1 unknown coefficients 
of Pn{x). We will choose the type II multiple orthogonal polynomials to be monic so 



Pfi{x)wi{x)x^ dx 
Pft{x)w2{x)x^ dx 



Pfl{x)Wr{x)x^ dx 
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that the remaining coefficients can be determined uniquely by the orthogonahty 
relations, provided the matrix of coefficients has full rank. 

In this paper the emphasis will be on type II multiple orthogonal polynomials. The unicity 
of multiple orthogonal polynomials can only be guaranteed under additional assumptions on 
the r weights. Two distinct cases for which the type II multiple orthogonal polynomials are 
given as follows. 

1. In an Angelesco system (Angelesco, 1918) the intervals Aj, on which the weights are 
supported, arc disjoint, i.e., Aj fl Aj = whenever i ^ j. Actually, it is sufficient that 

o 

the open intervals Ai are disjoint, so that the closed intervals Aj are allowed to touch. 

Theorem 1 In an Angelesco system the type II multiple orthogonal polynomial Pn{x) 
factors into r polynomials 0^=1 i^) , where each Qn. has exactly Uj zeros on A^ . 

Proof: Suppose Ph{x) has rrij < Uj sign changes on Aj at the points Xi, . . . , x^n.. Let 
Qm {x) — [x — Xi) ■ ■ - {x — Xm ), then Pfi{x)Qm {x) does not change sign on A^, and 
hence 

/ Pn{x)Qmj {x)Wj{x) dx ^ 0. 

But this is in contradiction with the orthogonality relation on A^. Hence Pfj(x) has 
at least nj zeros on Aj. Now all the intervals Aj {j — 1,2, ... , r) are disjoint, hence 
this gives at least \n\ zeros of Pfi{x) on the real line. The degree of this polynomials is 
precisely \n\, so there are exactly Uj zeros on each interval Aj. g 

2. For aji AT system all the weights are supported on the same interval A, but we require 
that the \n\ functions 

Wi{x),XWi{x), ... , x"'^~^Wi{x),W2{x),XW2{x), . . . , x"''^~^W2{x) , 

. . . , Wr{x),XWr{x), . . . , x'^''~^Wr{x) 

form a Chebyshev system on A for each multi-index n. This means that every hnear 
combination 

/■ 

^Qnj-l(x)Wj{x), 

with Qjij-i a polynomial of degree at most nj — 1, has at most — 1 zeros on A. 

Theorem 2 In an AT system the type II multiple orthogonal polynomial Ph{x) has 
exactly \n\ zeros on A. For the type I vector of multiple orthogonal polynomials, the 
linear combination ^j=i ^nj(a;)wj(a;) has exactly |n| — 1 zeros on A. 
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Proof: Suppose Pft{x) has m < \n\ sign changes on A at the points Xi, . . . .x^- Take 
a multi-index fh = {mi, 1712, ■ ■ ■ ,mr) such that rrii < rii for every i and rrij < rij for 
some j and construct the function 

r 

Qi^) = y^,Qiix)wi{x), 

where each Qi is a polynomial of degree mj — 1 whenever i ^ j, and Qj is a polynomial 
of degree m.j, satisfying the interpolation conditions 

Q{xk) = 0, k = 1,2, . . . ,m, 

and Q{xo) — 1 for an additional point xq G A. This interpolation problem has a 
unique solution since we are dealing with a Chebyshev system. The function Q has 
already m zeros, and since we are in a Chebyshev system, it can have no additional 
sign changes. Furthermore, the function does not vanish identically since Q{xo) = 1. 
Obviously Pfi{x)Q{x) does not change sign on A, so that 

/ Pri{x)Q{x)dx^^, 

but this is in contrast with the orthogonality relations for the type II multiple orthog- 
onal polynomial. Hence Pn{x) has exactly \n\ zeros on A. 

The proof for the type I multiple orthogonal polynomials is similar. First of all, since 
we are dealing with an AT system, the function 

r 

has at most \n\ — 1 zeros on A. Suppose it has m < \n\ — 1 sign changes at the 
points Xi, X2, . . . , Xm, then we use the polynomial Qm{x) = {x — Xi) ■ ■ ■ {x — Xm) so 
that A{x)Q{x) does not change sign on A, and 

/ A{x)Q{x)dxy^O, 

which is in contradiction with the orthogonality of the type I multiple orthogonal 
polynomial. Hence A{x) has exactly \n\ — 1 zeros on A. g 



Orthogonal polynomials on the real line always satisfy a three-term recurrence relation. 
There are also finite order recurrences for multiple orthogonal polynomials, and there are 
quite a few of recurrence relations possible since we are dealing with multi-indices. There is 
an interesting recurrence relation of order r+1 for the type II multiple orthogonal polynomials 
with nearly diagonal multi-indices. Let n G N and write it as n = kr + j, with < j < r. 
The nearly diagonal multi-index s{n) corresponding to n is then given by 

s{n) = i k + l,k + 'L... .k + LLk,.^. ,k) . 

j times r—j times 
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If we denote the corresponding multiple orthogonal polynomials by 



then the following recurrence relation holds: 

r 

^Pn{^) = Pn+l{^) + 0'n,jPn-j (^) ; 

j=0 

with initial conditions Po{x) — 1, Pj{x) — for j — —1, —2, .. . ,—r. The matrix 

/ ^0,0 1 \ 
ai,i ai,o 1 

^2,2 ^2,1 02,0 1 



(2.19) 



V 



Or 



n,r ^n,r— 1 



has eigenvalues at the zeros of P„+i (x) , so that in the case of Angelesco systems or AT systems 
we are dealing with non-symmetric matrices with real eigenvalues. The infinite matrix will 
act as an operator on i"^, but this operator is never sclf-adjoint and furthermore has not a 
simple spectrum, as is the case for ordinary orthogonal polynomials. Now there will be a set 
of r cyclic vectors and the spectral theory of this operator becomes more complicated (and 
more interesting). There are many open problems concerning this non-symmetric operator. 



3 Some very classical multiple orthogonal polynomials 

We will now describe seven families of multiple orthogonal polynomials which have the same 
flavor as the very classical orthogonal polynomials of Jacobi, Laguerre, and Hermite. They 
certainly deserve to be called classical since they have a Rodrigues formula and there is a first 
order differential operator which, when applied to these classical multiple orthogonal poly- 
nomials, gives another set of multiple orthogonal polynomials. However, these are certainly 
not the only families of multiple orthogonal polynomials (see Section 4.1). The first four 
families are AT systems which are connected by limit passages, the last three families are 
Angelesco systems which are also connected by limit passages. All these families have been 
introduced in the literature before. We will list some of their properties and give explicit 
formulas, most of which have not appeared earlier. 



7 



AT systems 



Angelesco systems 



Jacobi-Pineiro 

p(«0,"l,"2)/ \ 



Jacobi- Angelesco 

p(a,/3,7) 
J ri,m 




multiple Laguerre I 

r ("i."2)/ N 
J-'n.m \ J 



multiple Laguerre II 

r K,ci,c2)/ X 



Jacobi-Laguerre 




multiple Hermite 



3.1 Jacobi-Pineiro polynomials 

The Jacobi-Pifieiro polynomials are multiple orthogonal polynomials associated with an AT 
system consisting of Jacobi weights on [0, 1] with different singularities at and the same 
singularity at 1. They were first studied by Pineiro when = 0. The general case 
appears in ||3^, p. 162]. Let ao > —1 and ai, . . . , be such that each ai > —1 and 
ai — aj ^ Z whenever i ^ j. The Jacobi-Pineiro polynomial p^""'"-* for the multi- index 
n = (ni,n2, . . . ,n.r) G N*" and a = (ai, . . . , a^) is the monic polynomial of degree \n\ = 
rii + n2 + ■ ■ ■ + rir that satisfies the orthogonality conditions 

[ P^"'''"^(x)x°i(l = 0, = 0,1,... ,ni - 1, (3.20) 

Jo 



Pi'''""\x)x''^{l-xf'x''dx = 0, A; = 0,1,... ,712-1, (3.21) 







P^""'"^(x)x°'-(1 -x)"«x'=dx = 0, A; = 0,1,... ,n^-l. (3.22) 







Since each weight Wi{x) = x"*(l — a;)"° satisfies a Pearson equation 

x(l — x)w'^{x) = [ai{l — x) — aQx]wi{x) 
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and the weights are related by 



Wi{x) = x"' °''Wj{x), 



one can use integration by parts on each of the r integrals ( p.2C )-( 3.22| ) to find the following 
raising operators: 



d_ 

dx 



x-^{l - xrP^r'''\^) ] = + «o + - a:)""-^P^:r^'"-^^-(x), (3.23) 



-1 o("0-l,<3-ej) , 



where ej is the jth standard unit vector. Repeatedly using this raising operator gives the 
Rodrigues formula 



-1)1"! J](|n|+«o + «, + l)n,Pi' 



("0,a), 



a;) = (1 - x)-"" Yl 



dx'^i 



(3.24) 



The product of the r differential operators x~°'^ D^'^x^^'^"'^ on the right hand side can be 
taken in any order since these operators are commuting. 

The Rodrigues formula allows us to obtain an explicit expression. For the case r = 2 we 
write 

(-l)"+™(n + m + ao + tti + l)n(n + m + ao + as + l)mPt^"""'\x) 



dx^ 



'\ — nf,ai-a2+n ^a2+m^'^ _ ^-^ao+n.+m 



The mth derivative can be worked out using the Rodrigues formula ( |1.8|) for Jacobi polyno- 
mials and gives 



'-lY{n + m + ao + ai + l)„Pi™'"^Ha^) 



'1 - x)-"«a;-"i^a;°i+"(l - a;)"°+"P^"o+"'"=^)(a;). 



Now use Leibniz' rule to work out the nth derivative: 

= (1 - x)-''°x-'' 



dx^ dx""-^ 



In order to work out the derivative involving the Jacobi polynomial, we will use the following 
lemma. 

Lemma 1 Let Pi"'^\x) be the nth degree monic Jacobi polynomial on [0, 1]. Then for a > 
and P > —1 



d 



- [(1 - x)"Pi"'^)(x)] = -{a + n)il- xr~'P!C~''^^'\x), 



(3.26) 



and 



dx^ 



[(1 - x)"Pi"''^)(x)] = (-!)'"(« + n - m + 1)„(1 - x)"-'"Pi™'^+™)(a;). (3.27) 
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Proof: First of all, observe that 



^ [(1 - x)"p(°'^)(x)] = (1 - xr-' {-aPt^^\x) + (1 - x)[Pt^^\x)r) , 

so that the right hand side is —{a + n){l — x)°'^^Qn{x), with Qn a monic polynomial of 
degree n. Integrating by parts gives 

-{a + n) [ {l-xT-^x'^+^+^Qn{x)dx 
Jo 

_ ^/3+fc+l(l_^)apK/3)(2;)|;^-(/5 + fc + l) / x'^^'' {1 - x)'' Pj{'^^\x) dx. 



Obviously, when a > and (3 > —1, then the integrated terms on the right hand side vanish. 
The integral on the right hand side vanishes for k = 0,1, . . . ,n — l because of orthogonality. 
Hence Qn is a monic polynomial which is orthogonal to all polynomials of degree less than n 
with respect to the weight x^~^^{l — x)°'~^, which proves ( 3.26 ). The more general expression 
(|3.27|) follows by applying (|3.26|) m times. g 

By using this lemma we arrive at 

^n + m + ao + a^ + l)„Pi;^^'°^'"^)(a;) 

k=o ^ / V ^ / 
For the Jacobi polynomial we have the expansion 



(a 



+ (3 + n + l)„Pi"'^)(^) = ^' E t f) ^""'(^ - 1)'' (3-28) 



which can easily be obtained from the Rodrigues formula ( |1.8| ) by using Leibniz' formula, so 
that we finally find 



+ m + ao + ai + l)n(n + m + ao + "2 + l)^^^"^'"''"'^^;) 

ai + n\ f oq + m + n\ f a2 + n + m — k\ f ao + k + m\ ^ 



[n 

n m 

n\m\ 

k=0 j=0 



E E ( ) ( j y - ft. j l^ao ft. -t- in J ^n+m-k-j^^_-^^k+j 



(3.29) 



We can explicitly find the first few coefficients of Pm,n°'^'°'^\x) from this expression. We 
introduce the notation 

n!m! 

K„ 



(n + m + ao + ai + !)„(« + m + ao + a2 + l)m 
ao + ai + 2?T, + m\ ^ f ao + a2 + 2m + n" ^ 



n \ m 



First let us check that the polynomial is indeed monic by working out the coefficient of x™"'""'. 
This is given by 



fc=0 3=0 



ai + n\ f ao + m + n\ f a2 + n + m — k\ f ao + k + m 
k J\ n-k J\ j )\m-j 
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The sum over j can be evaluated using the Chu-Vandermonde identity 



E 

i=o 



a2 + n + m — k\ /"ao + k + m\ /ao + a2 + n + 2m 
3 )\rn-j) \ m 



which is independent of k. The remaining sum over k can also be evaluated and gives 

n 

E 



fc=0 



cti + n\ /ao + m + n\ / cto + cti + m + 2n 
k I \ n — k I \ n 



and the double sum is therefore equal to K^^^, showing that this polynomial is indeed monic. 
Now let us write 

p{ao,ai,a2) ( „\ _ „m+n , a n+m-1 , rj n+m-2 , n+m-3 _i_ . . . 

n,m \ ) — ~ ^n,m'*' ~ ^n,m-^ ~ ^n,m-^ ~ 

The coefficient An^^ of x"^~^"'~^ is given by 



n m / 
fc=0 j=0 ^ 



tti + /«o + + n\ /q;2 + ri + m — A;\ /ao + A; + m 
/c /\ n — k l\ i / \ m — i 



This double sum can again be evaluated using Chu-Vandermonde and gives 

n{ai + n){aQ + q;2 + n + m) + m{a2 + n + m)(ao + ai + 2n + m) 



(cco + CKi + 2n + m)(Q;o + q;2 + n + 2m) 



Similarly we can compute the coefficient Bn^rn of x""'""^"^ and the coefficient Cn,m of 

but the computation is rather lengthy. Once these coefficients have been determined, one 

can compute the coefficients in the recurrence relation 

xPn{x) = Pn+l{x) + hnPn{x) + C„P„_i(a;) + dnPn-2{x) , 

where 

P2n{x) = P^r'"^^\x), P2nM^) = 

Indeed, by comparing coefficients we have 

b2n — An^n ~ An+l^m ^2n+l = ^n+l,n ~ ^n+l,n+l ) (3.30) 

which gives 

b2n = [36n^ + (48q;o + 28q;i + 20q;2 + 38)n^ 

+ (21q;q + Sal + + ^Oagai + lSaoa2 + 15aia2 + S9ao + 19ai + 19a2 + 9)n 

+ (3q;q + 10q;qQ;i + 4q;qQ;2 + daoa^ + 2aoal + llaoaia2 + ha\a2 + Saia^ 

+ 12q;q + Sa^ + ?>a^ + 13q;oQ;i + 13q;oQ;2 + 8q;iq;2 + 6q;o + 3q;i + 30:2)71 

+ ccq + aoQii + q;2Q;^ + 2a20!^aQ + 2q;qQ;i + ctj^ao + ccgCCo + q;2Q;i + a^ai 

+ agQ;]^ + a^a^ai + 0:20;]^ + 2Q;2ttoQ;i + 3q;2Q;iq;o + 2q!2«o + ot\a2 + ckq + Q;oQ;2j 

X (3n + ao + Q;2)~^(3n + ao + Q;i)~^(3n + ckq + q;2 + l)~^(3n + ckq + ai + 2)"''", 
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and 

Wi = [sen'' + (48ao + 20ai + 28q;2 + 106)n^ 

+ {21aQ + Aaf + Sa^ + ISaoai + 30aoa2 + 15aia2 + lOSao + 41q;i + 65a2 + lll)n^ 

/322 22 2 2 

+ (Sttg + 4q;qQ;i + lQaQa2 + 2q;oQ;i + 6aoa2 + 11q;oQ;iq;2 + 3q;;^q;2 + 5aia2 

+ SOckq + 5al + ISaj + 23q;o«i + 47q;oQ;2 + 22q;iq;2 + 72ao + 25ai + A9a2 + 48)n 

+ 18aoQ;2 + 8q;2Q;o + 4q;i + Aalai + 8q;iQ;2 + 2q;q + ba^ao + 8Q;2aiao + 12a2 
+ 7 + 15q;o + C(2'^i + lOag + Qaoai + 2Q!2tti + 2q;qQ;i + a^ao + + a2aQ 

OOO O O O 1 

+ Q;2«o + + Ot2Ci\(y.Q + 2ci2(y.i^Cl\ + 2q;2«oChJ 

X (3n + cto + Q!2 + l)~^(3n + cto + cti + 2)~^(3n + + ct2 + 3)~^(3n + ao + cti + 3)""^. 
For the recurrence coefficient c„ we have the formulas 

C2n = -Bn,n ~ -Sn+l,n ~ ^2n^n,n) C2n+1 = -Bn+l,n ~ -Bn+l,n+l ~ ^2n+l^n+l,n) (3.31) 

which, after some computation (and using Maple V), gives 

C2n = n{2n + ao) (2n + ao + «i) (2n + ao + ci!2) 
X [54n'' + (63ao + 45ai + 45a2)n^ 

+ (24ao + Sal + 8^2 + 42aotti + 42aott2 + 44aia2 - 8)n^ 

,0000 Q 00 

+ (3q;q + + a2 + 12aQai + 12q;qQ;2 + 3q;oQ;i + 3q;oQ;2 + 33q;oQ;iQ;2 + 8q;iQ;2 
+ 8aia2 — 3^0 — 4q;i — Aa2)n 

+ aQfti + aQQ.2 + 6aQaia2 + 0^10^2 + cnck2 + 3aoaiCt2 + 3ao<^i<^2 ~ ctofti — aoC(2 — 2ci\OL2\ 
X (3n + ao + «! + l)^^(3n + ao + <^2 + l)^^(3n + cto + Q;i)^^(3n + cto + ^2)"^ 
(3n + tto + t^i ~ l)^^(3n + cto + 0:2 — 1)^^ 

and 

C2n+i = (2n + ao + 1) (2n + ao + cti + 1) (2n + + q;2 + 1) 
X [54n^ + (63ao + 45ai + 45^2 + 135)n^ 

+ (24^0 + Sa\ + 8^2 + 42aoai + 42aoa2 + 44aia2 - +126ao + 76ai + 104q;2 + 120)n^ 

,0000 Q 00 o 

+ (3aQ + + + 12q;qQ;i + 12q;qQ;2 + 3q;oQ;]^ + 3q;o<^2 + 33ao'^iCt2 + 8a|Q;2 
+ 8aia2 + 36«Q + So;^ + 19q;2 + +54Q;oQii + 72Q;oQi2 + 66q;iQ;2 + 87q;o + 39q;i 
+ 81^2 + 45)n2 

+ (ag*^! + ao'^2 + 6aQaia2 + ttia2 + ci\0i2 + 3ao«i'^2 + 3ao<^iCi2 + 3tto + 2«2 



+ 12Q;Qai + \2aQOL2 + 6aoa2 + 33aoaia2 + hoL^OL2 + WdxQ.^ + 18q;o + 20q;oQ;i 
+ 38aoCt2 + 14«2 + 26aia2 + 24ao + 6ai + 24q;2 + 6)n 

+ ttQai + 3q;qQ;iQ;2 + 3q;oQ;iQ;2 + ttiQ;2 + + 0:2 + 3Q;Qai + 3q;qQ;2 + 'ociQOL\OL2 
+ 3q;o«2 + 3aiQ;2 + 3q;q + 3^2 + 2«oCH + 6q;oQ;2 + 2q;iq;2 + 2q;o + 20:2] 
X (3n + tto + tti + 3)~^(3n + cto + ^2 + 2)~"'(3n + ao + cti + 2)~^(3n + ao + 0^2 + 
(3n + ao + ai + l)"-^(3n + ao + a2)~^ 

Finally, for dn we have 

^2n — Cn,n Cn+l,n ^2n-^n,n C2n^n,n— 1; 
C^2n+1 — Cn+l,n ~ Cn+l,n+l ~ ^2n+l-Sn+l,n ~ C2n+l-^n,m (3.32) 
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giving 

d2n 



n{2n + aQ){2n + ao — l){2n + ao + ai)(2n + ao + ai — 1) 

(2n + tto + 0:2) (2n + + «2 -!)('"' + + ai - 02) 

(3?2 + 1 + tto + ai)"^(3?2 + tto + ai)^^(3n + ao + ^2) ""^(S^^ — 1 + + 

(3n — 1 + tto + a;2)~"'^(3r7, — 2 + ao + ai)~^{?>n — 2 + ao + 0.2)'^ 



and 

d2n+l 



n{2n + 1 + aQ){2n + aQ){2n + ao + ai)(2n + 1 + ao + O-i) 
(2n + 1 + ao + a2)(2n + ao + a2)('"' + OL2){n + ^2 - "i) 

(3n + 2 + ao + ai)~^{?)n + 2 + ao + a;2)^"'^(3r2 + 1 + ao + Q;i)^^(3n + 1 + ao + «2) 
(3?T, + ao + ai)~^(3n + ao + a;2)~^(3n — 1 + + tt2)~^- 

These formulas are rather lengthy, but explicit knowledge of them will be useful in what 
follows. Observe that for large n we have 



-2 



lim hn 

n— >oo 

lim Cn 

n— >oo 

lim (in 



4 _ 
9 ~ 
16 
243 
64 
19683 



4 

27 



4 

27 



3.2 Multiple Laguerre polynomials (first kind) 

In the same spirit as for the Jacobi-Pifieiro polynomials, we can consider two different families 
of multiple Laguerre polynomials. The multiple Laguerre polynomials of the first kind L'i{x) 
are orthogonal on [0, 00) with respect to the r weights Wj{x) = x°'^e~^, where aj > —1 for 
j = 1,2,... , r. So these weights have the same exponential decrease at 00 but have different 
singularities at 0. Again we assume — a-,- ^ Z in order to have an AT system. These 
polynomials were first considered by Sorokin |?T[]. The raising operators are given by 



)) 



-x^'^-^e- 



J 



1, 



and a repeated application of these operators gives the Rodrigues formula 



-i)I"IlS( 



n 



X 



^3 



dx^' 



(3.33) 



(3.34) 



When r = 2 one can use this Rodrigues formula to obtain an explicit expression for these 
multiple Laguerre polynomials, from which one can compute the recurrence coefficients in 



xPn{x) = Pn+l{x) + bnPn{x) + C„P„_i(a;) + dnPn-2{x 

where P2n{x) = L'n^n°''^\x) and P2n+l{x) = L^n+l,n 
Jacobi-Pifieiro polynomials, it is much easier to use the limit relation 



L)^li"^'(x). But having done all that work for 



(3.35) 
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The recurrence coefficients can then be found in terms of the following limits of the corre- 
sponding recurrence coefficients of Jacobi-Pineiro polynomials: 



bn 



lim 6(fO,ai,«2) 



ao, 



lim ci°0'"^'"^)a^, 
lim ^"'''"''"'^"o. 



givmg 



b2n = 3ri + «! + 1, 

&2n+l = 3ra + 02 + 2, 

C2„ = n(3n + ai + a2), 

C2n+i = Sn^ + (ai + a;2 + 3)n + ai + 1, 

= n{n + ai){n + ai - a2), 

c?2n+i = n{n + a2){n + a2 - ai). 



Observe that for large n we have 



bn 

lim — 

n— >cxa ?2 



1 . 

n^oD fl^ 
dn 

lim — 

n— >oo 77,'' 



2 V2 
4 V2 



1 



3.3 Multiple Laguerre polynomials (second kind) 

Another family of multiple Laguerre polynomials is given by the weights Wj{x) = x"'"e~'^^^ 
on [0, oo), with Cj > and Cj ^ Cj for i ^ j. So now the weights have the same singularity 
at the origin but different exponential rates at infinity. These multiple Laguerre polynomials 
of the second kind L^°'''^(x) appear already in p. 160]. The raising operators are 



dx 



("0,c) 



(3.36) 



and a repeated application of these operators gives the Rodrigues formula 



r("o,c) 



XI = X 



n 



dx"'^ 



X 



\n\+ao 



(3.37) 



These polynomials are also a limit case of the Jacobi-Pineiro polynomials. For the case r = 2 
we have 



LlT"''\x) = lim (-a)"+™p("r^"'^^°)(l - x/a). 



(3.38) 
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The recurrence coefficients can be obtained from the corresponding recurrence coefficients 
of Jacobi-Pineiro polynomials by 

bn = lim (1 - 6(,"°''^i"'^2"))q;, 

a— >oo 

Cn = lim C^"0''^i'*''^2a)Q,2^ 

a— >oo 

dn = lim 



givmg 



->2n+l 



C2n 



n{ci + 3C2) + C2 + q;oC2 

C1C2 

n(3ci + C2) + 2ci + C2 + ttoCi 
C1C2 

n(2n + q;o)(ci + cl) 

2 2 ' 
C1C2 



_ 2v?{c\ + c|) + n[cl + 3c^ + Q;o(cf + c|)] + c| + ctoci 

C2n+1 — 2"^ ' 

_ n(2n + Q!o)(2n + ao - l)(c2 - Cl) 

W2n — 3 1 



qc2 

n(2n + Q;o)(2n + ao + l)(ci - C2) 



C^2n+1 — 3 



Observe that for large n we have 



lim 

n— >oo n 



I Cl + 3C2 _ 

L it n = U (mod 2), 



2CiC2 
3Ci + C2 

2CiC2 



if n = 1 (mod 2), 



Cn ^1 + c| 



lim 



I '^^ o '^^ if n = (mod 2), 
3 if n = 1 (mod 2). 



1 "^2 



2ciC2 

3.4 Multiple Hermit e polynomials 

Finally we can consider the weights Wj{x) — e^^^+'^j^ on (—00, 00), for = 1, 2, . . . , r and 
Cj different real numbers. The multiple Hermite polynomials H^{x) once more have raising 
operators and a Rodrigues formula, and they are also limiting cases of the Jacobi-Pihciro 
polynomials, but also of the multiple Laguerre polynomials of the second kind. For r = 2 
this is 

H^'l'^'\x) = lim (2^)"+"^p(«."+'^iV^.°+-2V^) ( ^ + ^ \ ^ (3 39) 
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so that the recurrence coefficients can be obtained from the Jacobi-Pineiro case by 

a— >oo 2 



c„ = hm 4c("."+^i^."+^2v^) 



a 



3 



This gives 



hn 


= Ci/2, 




'2n+l 


= C2/2, 




Cn 


= ^/2, 




d2n 


= ^(ci - 


C2)/4, 


''2n+l 


= n{c2 - 


ci)/4. 



Alternatively, we can use the limit transition from the multiple Laguerre polynomials of the 
first kind: 



^^rn^) = lim a"+-LSr ^^"/''""'"^"/'^v^a: + a). (3.40) 

' a—^oo 

The recurrence coefficients are then also given in terms of the following limits of the recur- 
rence coefficients of the multiple Laguerre polynomials of the first kind 



I. ^■ /,("+ci-\/a72,a+C2A/a72) x , 

On = hm (o„ ^ ^ -a)/v2a 



(ci+ciA/a72,a+C2A/a72) ,/„ ^ 

c„ = hm Cn /(2q;), 

(a+ciya72,a+C2-v/a72) ,, /7^x3 

dn = lim dn /(v2«) , 

which leads to the same result. Observe that for large n we have 

K 



lim — ^ = 



lim 

n^oo n 2 

3.5 Jacobi-Angelesco polynomials 

The following system is probably the first that was investigated in detail (||2^, ||2^)- It is an 
Angelesco system with weights Wi{x) = \h{x)\ on [a, 0] (with a < 0) and W2{x) = \h{x)\ on 
[0, 1], where h{x) = (x — a)"x^(l — xy and a, 7 > —1. Hence the same weight is used for 
both weights wi and W2 but on two touching intervals. The Jacobi-Angelesco polynomials 
PrCm'"'^ {x; a) therefore satisfy the orthogonahty relations 

P!Cm'^\x;a){x-ay\xf{l-xyx^dx = 0, A; = 0, 1, 2, . . . , n - 1, (3.41) 

1 

P^"^^'^)(x;a)(x-a)"x'^(l-x)Vda; = 0, = 0, 1, 2, . . . , m - 1. (3.42) 
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The function h{x) satisfies a Pearson equation 

(x — a)x(l — x)h\x) = [ax{l — x) + [3{x — a)(l — x) — 7(x — a)x]/i(x), 

where {x — a)x{l — x) is now a polynomial of degree 3. Using this relation, we can integrate 
the orthogonality relations by part to see that 



d 

dx 



[{x-arx^{l-xr>Pi(:^^^\x-a)\ 



-{a + P + -i + n + m){x- ay-^x^-'il - x^-' Pil-lf;,''^~'\x; a), (3.43) 



which raises both indices of the multi-index {n,m). Repeated use of this raising operation 
gives the Rodrigues formula 

d"" 



dx^ 



= + (3 + ^ + k + 2m + l)m{x - a^x^il - Pj^f,]l{x; a). (3.44) 

For k = and m = n, this then gives 

[{x - a)°+"x^+"(l - x)^+"] 

= (-l)"(a + /5 + 7 + 2n + l)„(x - a)"x''(l - x)^ Pi'^/'^Xx; a). (3.45) 

Use Leibniz' formula to find 



(-l)"(a + (3 + ^ + 2n + l)„(x - a)"a;'^(l - x)^P^°/'^)(x; a) 

fc=0 



A; / V dx'' 



dx^ ^ 



Now use the Rodrigues formula for the Jacobi polynomials ( |1.8| ) to find 



E(-ir'( 



/5 + 7 + 2n\ fa + n 



k 



n — k 



(x - a)'=x"-'(l - x)"~'=P,(^+"-''^+"-')(x). 



Use of the expansion (|3.28| ) for the Jacobi polynomial gives 



a 



+ ^ + ^ + 2^Vi-^-)(x;a) 



n 

n k 



j n,n 



fc=0 j=0 

n n—k 



n — kj \ j J \k — j 

7 + 
n — k — j 



E g (" I ") f ; " I (... ' t " ,■ ) - 1)'*'- 



(3.46) 
(3.47) 
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where the last equation follows by the change of variable k t-^ n — k. If we write this in 
terms of Pochhammer symbols, then 



a + P + '-y + 3n 



n 



(7 + l)n {-n)k+j{-a-n)k{-(3-n)j _ .n-k( _ i)k+j n-j 

nr \nfl + ^\ t:^ f n X - 1 X - 1 

x {x-aj { Fi -n, -a - n, -/5 - n, 7 + 1 



T7j / \ 3^ (Xi 



Fi(a, b,b',c;x,y) = 



where 

is the first of Appell's hypergeometric functions of two variables. 
For the polynomial P^"'f {x; a) we have the Rodrigues formula 

Jn p 



dx^ 

= (-!)"(« + /? + 7 + 2n + 2)„(x - a)°x^(l - x)^Pi+f a), (3.48) 

where p^('^+"'^+"-''T+") j^^;- _ 2; _ X^'^''^'' is the monic orthogonal polynomial of first degree 
for the weight (x — a)"'^"|a;|^^"(l — x)"'^^ on [a, 0]. If we write down the orthogonality of 
this polynomial to the constant function, 



(x - X("''^'^))(x - a)"+"|x|^+"(l - xy^"" dx = 0, 



then we see that 

rO 



£ x{x - a)"+"|a;|^+"(l - x)"'+'' dx 



{x - a)°+"|x|^+"(l - dx 
A standard saddle point method gives the asymptotic behavior 

lim X^"'^'^) = xi, (3.49) 

n— >oo 

where Xi is the zero of g'{x) in [a, 0], where g{x) = {x — a)x{l — x). . Combining the 
Rodrigues equation in ( ^.481 ) with the Rodrigues equation ( |3.45| ) shows that 

Pif^::\x; a) = xPi^r'^^\x- «) - ^^'^^^^f^^^f^^ir/^Ha:; a). (3.50) 
In order to compute the coefficients of the recurrence relation 

xPn{x) = Pn+lix) + hnPn{x) + C„P„-l(x) + rf„P„_2(x), 

where 

P2n{x) = p!Ci'''\x; a), P2„-,i(x) = Pi^f a). 
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we will compute the first few coefficients of the polynomials 

p(">/3.7)('^. — rf-m+n _|_ A n+m-l , n m+n-2 , ^ n+m-3 _i_ . . . 

First we take n = m. In order to check that our polynomial is monic, we see from ( p.46|) 
that the leading coefficient is given by 

a + P + •y + 3n\ f a + n\ f P + n\ f •y + n 



n I ^ ^ \n — k J \ i J \ k — j 



Chu-Vandermonde gives 

k 

E 

3=0 



[3 + n\h + n\ _ f(3 + y + 2n 



and also 



E 



a + n\//? + 7 + 2n\ _ / a + /3 + y + 3n 
n — k I \ k I V n 



k=0 

J2n~\ 



SO that the leading coefficient is indeed 1. The coefficient An.n of ^ is equal to 



— 1 71 k 



„ , EE r;- r.-x--— 



a + /5 + 7 + 'in\ ( a \ n \ ( P \ n \ ( y \ n 

k=o j=o 

Working out this double sum gives 



^(aA7) = -n[a + l3 + 2n + a{P + y + 2n)] ^ ^^^^^^ 
a + P + y + 3n 



For 

+1 n coefficient v4„+i „ of x^" can be obtained from ( p.50| ) 

= 4^^^'^^ - + ^ + ^ + f + (3.52) 

The coefficient 6„ in the recurrence relation can now be found from (p.30|) 

^ ^ n[n + 7 + a(n + a)] ^^(^^^^^) 2n + a + P + y + 1 

(a + /5 + 7 + 3n)(a + /? + 7 + 3n + l) " Sn + a + P + y + V 

b2n+i = (^5n2 + (4a + 4/? + 37 + 7)n + (a + /5 + 7 + l)(a + /? + 2) 

+ a[5n^ + (3a + 4/5 + 47 + 7)n + (a + /? + 7 + + 7 + 2)]) 

X {a + P + y + 3n + iy^{a + p + y + 3n + 3)"^ 
2n + a + P + y + l 
3n + a + P + 'j + l' 
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The coefficient B^^n of x^" ^ in P^yf''^\x] a) is given by 
^(a,/3,7) ^ an(Q + /3 + 7 + 2n)(/? + n) 

(a + /? + 7 + 3n)(a + /5 + 7 + 3n - 1) 
n(n — 1) 



+ 



2(a + /3 + 7 + 3n)(a + /3 + 7 + 3n - 1) 
X [(a + /5 + 2n)(a + /3 + 2n - 1) + 2a(a + /3 + 2n)(/5 + 7 + 2n) 
+ a\(3 + 7 + 2n)((/5 + 7 + 2n - 1)] , 
and from (|3.50| ) we also find 

T^(a,/3,7) ^ n(a,/3+l,7) _ y("./^.7) 4(q./3.7) « + /^ + 7 + 2n+l 
n+l,n n,n n n,n i o i i o i i ' 

Q! + p + 7 + 3n + l 



Using ( p.31|) then gives 

n{a + /3 + 7 + 2n) 



C2ri 



(a + /? + 7 + 3n-l)(a + /5 + 7 + 3n)2(a + + ^ + 3,2 _ 1) 
X ((a + /? + 2n)(7 + n) - 2a{a + n)(7 + n) + a^(/3 + 7 + 2n)(a + n) ) , 



and 

C2n+1 



a + /? + 7 + 2ri + l 



(a + /5 + 7 + 3n + 3)(a + /5 + 7 + 3n + 2)(a + /3 + 7 + 3n + l)2(a + + 7 + 3n) 
X (^n(n + 7)(a + /? + 2?2 + l)(a + /? + 7 + 3n + 3) 
- a[24n^ + (29a + 41/5 + 297 + 48)^^ 

+ (lOa^ + 39a/3 + 26^7 + 29/3^ + 39/^7 + 107^ + 44a + 62/5 + 447 + SO)^^ 
+ (a^ + lW(3 + 5a27 + 19a/?2 + 24a/?7 + 5a72 + 9/?^ + 19/5^7 + ll(3-i'^ + 7^ 
+ lla^ + 39a/3 + 28a7 + 28/3^ + 39/^7 + 117^ + 19a + 25/3 + I97 + 6)n 
+ (a + /3 + 7)(a + /3 + 7 + l)(a + /3 + 7 + 2)(/3 + 1)] 

+ a?n{n + a)(/3 + 7 + 2n + l)(a + /? + 7 + 3n + 3)' 



, a + l3 + l + 2n + l ^(^^^^^^ 

(a + /? + 7 + 3n + 3)(a + /? + 7 + 3n + l)2(a + /3 + 7 + 3n) " 

X (^12^3 + (16a + 16/3 + IO7 + 18)71^ 

+ [{a + l3 + 7) (7a + 7/3 + 27) + 16a + 16/3 + 107]n 
+ (a + /3 + 7)2(a + /3) + (a + /3 + 7)(3a + 3/3 + 27 + 2) 
+ a [I2n^ + (10a + 16/3 + I67 + 18)n2 
+ [(a + /3 + 7) (2a + 7/3 + 77) + 10a + 16/3 + 167]n 
+ (a + /3 + 7)2(/3 + 7) + (a + /3 + 7)(2a + 3/3 + 37 + 2)]) 
_ (a+^+7 + 2n+2)^^^(„,^,^)^2_ 



(a + /3 + 7 + 3n + 1)2 

The coefficient C„,„ of in Pi^'^'^''(x; a) can be computed in a similar way, and the 

coefficient Cn+i,n of x^""^ p^"'^^)(x; a) is given by 

^(a,/3,7) ^ ^(a,/3+l,7) _ v(a,P,-y) p(a,/3,7) Q + /^ + 7 + 2^2 + 1 
W+l,n W,n n n,n ^ + ^ + ^ + 3^ + ^ • 
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A lengthy but straightforward calculation, using (|3.32|) , then gives 

^ _ -an{n + (3){a + 13 + -f + 2n){a + (3 + 'y + 2n - l)[n + -f + a{n + a)] 

^'^ ~ (a + /3 + 7 + 3n-2)(a + /5 + 7 + 3n-l)(a + /3 + 7 + 3ny{a + /5 + 7 + 3n + 1) 

n(a + /3 + 7 + 2n)(a + /3 + 7 + 2n - l)X^"f' 

(a + /5 + 7 + 3n-2)(a + /? + 7 + 3?2-l)(a + /? + 7 + 3?2)2(a + /5 + 7 + 3n + 1) 
X [(ri + 7)(a + /? + 2n) - 2a(n + 7)(?2 + a) + a^(n + «)(/? + 7 + 2n)] , 



and 

d2n+l 



n(a + /5 + 7 + 2n + l)(a + /3 + 7 + 2n) 



(a + /3 + 7 + 3n + 2)(a + /3 + 7 + 3n + l)2(a + /3 + 7 + 3n)^{a + /3 + 7 + 3n - 1) 
X (^(n + 7)(a + /3 + 2n)(a + /5 + 2ra + l) 

- a{n + a){n + -f){2a + 2/3 - 7 + 3n + 1) 

- a^(n + a)(n + 7)(-a + 2/5 + 27 + 3n + 1) 

+ a3(n + a)(/? + 7 + 2n)(/? + 7 + 2n + l)) 

n(Q + ;g + 7 + 2n + + ;g + 7 + 2n)xi"'^'^^ 

(a + /5 + 7 + 3n + l)2(a + /3 + 7 + 3n)2(a + /? + ^ + 3n - 1) 
X [{n + -f){a + (3 + 2n) - 2a{n + + 7) + a^(ra + a)(/5 + 7 + 2n)] . 

The asymptotic behavior of these recurrence coefficients can easily be found using ( ^.49| ), 
giving 

, a + 1 2x1 , 5(a + l) 2xi 

lim = — h ^— , lim 62n+i - 



n^oo 



3; "^;t-ri ri 

n^oo y 



lim C2n = -^(a^ - a + 1), lim C2„+i = -^x^ + J-Xi + ^(4a^ - a + 4), 

4 

lim (i2„ = — [2(a^ - a + l)xi - a(a + 1)], 

n— >oo ^4^5 

lim d2n+i = T^i^a^ - Sa' - 3a + 4) - ^(a^ - a + 1), 

where xi is the zero of g'{x) in [a, 0] and g{x) = {x — a)x{x — 1). These formulas can be 
made more symmetric by also using the zero X2 of g'{x) in [0, 1] and using the fact that 
xi + X2 = 2(a + l)/3: 

, a + 1 2xi , a + 1 2x2 

hm b2n = + 1™ f^2n+l = 5-, 

lim c„ = -^(a^ - a + 1), 

4 4 
lim c?2„ = -— 5'(xi), lim = -T;zg{x2). 
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3.6 Jacobi-Laguerre polynomials 

When we consider the weights wi{x) = (x — a)°\x\^e~^ on [a,0], with a < 0, and W2{x) = 
(x — a)"|x|^e~^ on [0, oo), then we are again using one weight but on two touching intervals, 
one of which is the finite interval [a, 0] (Jacobi part), the other the infinite interval [0, oo) 



(Laguerre part). This system was considered by Sorokin [38|. The corresponding Jacobi- 
Laguerre polynomials L^njn (x; a) satisfy the orthogonality relations 



^~^Tl 7YI 1 Ci) ^X Ci) X C X diX 



0, 



A; = 0,1, 



= 0,1 



, . . . , 



m — 1. 



The raising operator is 



^ [(x + a)"x^e-^L(:^;f)(x; a)] = -(x - a)"-^x^-^e-^xti^;i^)(x; a), 
from which the Rodrigues formula follows: 



[X 



- a)"+"'x'^+'"e-^L 



{a+m,l3+m) 
k,0 ' 



x; a] 



-irix - a)°x^e-"L^"''^^ 



(3.53) 



(3.54) 



From this Rodrigues formula we can proceed as before to find an expression for the polyno- 
mials, but it is more convenient to view these Jacobi-Laguerre polynomials as a limit case 
of the Jacobi- Angelesco polynomials 



Ll^,i\x;a) = hm 7"+-Pifi'^)(x/7; a/7), 



(3.55) 



so that ( |3.47D gives 



n n—k 



k=0 j=0 



a + n\f(3 + n\ (-l)'^+-'(x - a)"-'=x"-^' 
k 



(3.56) 



J / {n-k-j)\ 
For the recurrence coefficients in 

xPn{x) = P„+l(x) + fe„P„(x) + C„P„_i(x) + rf„P„_2(x), 

where P2n{x) = L^n,n\x] a) and P2n+i{'x) = L''^^^\{x; a) we have in terms of the correspond- 
ing recurrence coefficients of the Jacobi-Angelesco polynomials 



hm 76!"'^'^) (a/7), 

y— »oo 

hm fci"'^'^na/7), 
hm 7'4"'^'^n«/7), 



7^00 



and 



lim 7X, 

7— »oo 



(a,/3,7)| 
n 



ah) 



/° x(x — a)""'""|x|^"'""e ^ dx 



/°(x - a)"+"|x|^+"e-^' cix 
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:= X^°''^\ 



This gives 

b2n = n + X^^\ 

62n+i = 3n + a + /5 + 2 + a-X("'/^), 

C2n = n{a + P + 2n), 

C2n+i = n(« + /5 + 2n + l)-a(n + /3 + l) + (« + /3 + 2n + 2 + a)X^^)-(X("'^))2, 

4n = -an(/5 + n) +n(a + /3 + 2n)xi"'f , 
d2n+i = + + 2n)(a + + 2n + 1) + a(n + a)] - n(a + /3 + 2n)X^°'^\ 

For large n we have X^'^^ = a/2 + o(l) so that 
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3.7 Laguerre-Hermite polynomials 

Another hmit case of the Jacobi-Angelesco polynomials are the multiple orthogonal polyno- 
mials Hn%{x) for which 

/ H^^l{x)\xfe-''\''dx = 0, A; = 0,1, ... ,n- 1, 

J — oo 

POO 

/ H^%{x)x^e-'''x'' dx = 0, = 0, 1, . . . , m - 1. 

Jo 

We call these Laguerre-Hermite polynomials because both weights are supported on semi- 
infinite intervals (Laguerre) with a common weight that resembles the Hermite weight. These 



polynomials were already considered (for general r) by Sorokin [4C]. The limit case is ob- 
tained by taking 

Hi%{x) = lim (v/^)"+'"p(f^^'")(a;/v/^; -1). (3.57) 

This allows us to obtain the raising operator, the Rodrigues formula, an explicit expression, 
and the recurrence coefficients by taking the appropriate limit passage in the formulas for 
the Jacobi-Angelesco polynomials. For the recurrence coefficients this gives 

6„ = lim v^6i"'^'")(a = -1), 

Cn = lim ac(,"'^'")(a = -1), 

dn = lim(v^)=^4"''^'")(a = -l). 



and 



hm v^X("'^-")(a = -1) = ^ := ), 

J — oo I I 
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from which we find 



For large n we have 



so that 





= 








= n/2, 
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-1/2 ifn = (mod 2), 
1/2 if n = 1 (mod 2), 



1/4, 



-1/8 ifn = (mod 2), 
1/8 ifn = l (mod 2). 



4 Open research problems 

In the previous sections we gave a short description of multiple orthogonal polynomials and 
a few examples. For a more detailed account of multiple orthogonal polynomials we refer to 
Aptekarev Q and Chapter 4 of the book of Nikishin and Sorokin [Q. Multiple orthogonal 
polynomials arise naturally in Hermite-Pade approximation of a system of (Markov) func- 
tions. For this kind of simultaneous rational approximation we refer to Mahler ||2^ and de 
Bruin 0, [|T^. Hermite-Pade approximation goes back to the nineteenth century, and many 
algebraic aspects have been investigated since then: existence and uniqueness, recurrences, 
normality of indices, etc. The more detailed analytic investigation of the zero distribution, 
the nth root asymptotics, and the strong asymptotics is more recent and mostly done by 
researchers from the schools around Nikishin |3^, [33] and Gonchar |T^, [0. See in par- 
ticular the work of Aptekarev 0], Kalyagin [^, Bustamante and Lopez [^, but also 
the work by Driver and Stahl [jl5l, [T^ and Nuttall First one needs to understand the 
analysis of ordinary orthogonal polynomials, and then one has a good basis for studying this 
extension, for which there are quite a few possibilities for research. 



4.1 Special functions 

The research of orthogonal polynomials as special functions has now led to a classification 
and arrangement of various important (basic hypergeometric) orthogonal polynomials. In 
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Section 3 we gave a few multiple orthogonal polynomials of the same flavor as the very 
classical orthogonal polynomials of Jacobi, Laguerre, and Hermite. Regarding these very 
classical multiple orthogonal polynomials, a few open problems arise: 

1. Are the polynomials given in Section 3 the only possible very classical multiple or- 
thogonal polynomials? The answer very likely is no. First one needs to make clear 
what the notion of classical multiple orthogonal polynomial means. A possible way is 
to start from a Pearson type equation for the weights. If one chooses one weight but 
restricted to disjoint intervals, as we did for the Jacobi- Angelesco, Jacobi-Laguerre, 
and Laguerre- Hermite polynomials, then Aptekarev, Marcellan and Rocha |^ used the 
Pearson equation for this weight as the starting point of their characterization. For 
several weights it is more natural to study a Pearson equation for the vector of weights 
{wi, W2, . . . , Wr). Douak and Maroni [0, fl^ have given a complete characterization 
of all type II multiple orthogonal polynomials for which the derivatives are again type 
II multiple orthogonal polynomials (Hahn's characterization for the Jacobi, Laguerre, 
and Hermite polynomials, and the Bessel polynomials if one allows moment function- 
al which are not positive definite). They call such polynomials classical (i-orthogonal 
polynomials, where d corresponds to our r, i.e., the number of weights (functionals) 
needed for the orthogonality. Douak and Maroni show that this class of multiple 
orthogonal polynomials is characterized by a Pearson equation of the form 

+ = 0, 

where w = {wi, . . . ,WrY is the vector of weights, and \1/ and $ are r x r matrix 
polynomials: 



/ 1 ■■ 
2 ■■ 





\tpix) Ci C2 



\ 





r — 1 

Cr-l J 

with iIj{x) a polynomial of degree one and Ci, . . . , c^-i constants, and 



(<pi,i{x) 0i,2(a;) 
02,1 (a;) 02,2 (a;) 



01, r(a;)\ 

02, r(a;) 



\0r,l(a;) 0r,2(a;) ■■■ (f)r,r{x) J 



where 4>ij{x) are polynomials of degree at most two. In fact only 0^,1 can have degree at 
most two and all other polynomials are constant or of degree one, depending on their 
position in the matrix $. Douak and Maroni actually investigate the more general 
case where orthogonality is given by r linear functionals, rather than by r positive 
measures. We believe that Hahn's characterization is not the appropriate property to 
define classical multiple orthogonal polynomials, but gives a more restricted class. None 
of the seven families, given in the present paper, belong to the class studied by Douak 
and Maroni, but their class certainly contains several interesting families of multiple 
orthogonal polynomials. In fact, the matrix Pearson equation could result from a 
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single weight (and its derivatives) satisfying a higher order differential equation with 
polynomial coefficients. As an example, one can have multiple orthogonal polynomials 
with weights Wi{x) = 2x"+''/'^K^{2^/^) and ^2(0;) = 2x"+'^''+^^/'^K^+i{2^/^) on [0,oo), 
where Ky{x) is a modified Bessel function and a > — 1, z/ > (see |^ and [p!2| ). 

The polynomials of Jacobi, Laguerre, and Hermite all satisfy a linear second order dif- 
ferential equation of Sturm-Liouville type. A possible way to extend this characterizing 
property is to look for multiple orthogonal polynomials satisfying a linear differential 
equation of order r + 1. Do the seven families in this paper have such a differential 
equation? If the answer is yes, then an explicit construction would be desirable. We 
only worked out in detail the case where r = 2, so the search is for a third order 
differential equation for all the polynomials considered in Section 3. Such a third order 



equation has been found for certain Jacobi- Angelesco systems in ||2^. For the Angele- 
sco systems in Section 3 this third order differential equation indeed exists and it was 
constructed in |^. The existence (and construction) is open for the AT systems. A 
deeper problem is to characterize all the multiple orthogonal polynomials satisfying a 
third order (order r + 1) differential equation, extending Bochner's result for ordinary 
orthogonal polynomials. Observe that we already know appropriate raising operators 
for the seven systems described in Section 3. If one can construct lowering operators 
as well, then a combination of the raising and lowering operators will give the differ- 
ential equation, which will immediately be in factored form. Just differentiating will 
usually not be sufficient (except for the class studied by Douak and Maroni): if we 
take P^„(a;), then this is a polynomial of degree n + m — 1, so one can write it as 
Pn-i,m{,x)+ lower order terms, but also as Pn,m-i{,x)+ lower order terms. So it is not 
clear which of the multi-indices has to be lowered. Furthermore, the lower order terms 
will not vanish in general since there usually are not enough orthogonality conditions 
to make them disappear. 

3. In the present paper we only considered the type II multiple orthogonal polynomi- 
als. Derive explicit expressions and relevant properties of the corresponding vector 
{An^m{.x), Bn^m{.x)) of type I multiple orthogonal polynomials. Type I and type II 
multiple orthogonal polynomials are connected by 



Pn.mix) = const. 



but from this it is not so easy to obtain the type I polynomials. 

4. So far we limited ourselves to the very classical orthogonal polynomials of Jacobi, 
Laguerre, and Hermite. Discrete orthogonal polynomials, such as those of Charlier, 
Kravchuk, Meixner, and Hahn, can also be considered and several kinds of discrete 
multiple orthogonal polynomials can be worked out. It would not be a good idea to 
do this case by case, since these polynomials are all connected by limit transitions, 
with the Hahn polynomials as the starting family. At a later stage, one could also 
consider multiple orthogonal polynomials on a quadratic lattice and on the general 
exponential lattice, leading to g-polynomials. Again, all these families are related, 
with the Askey- Wilson polynomials as the family from which all others can be ob- 
tained by limit transitions. Do these polynomials have a representation as a (basic) 
hj^ergeometric function? Recall that we needed an Appell hypergeometric function of 
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two variables for the Jacobi-Angelesco polynomials, so that one may need to consider 
(basic) hypergeometric functions of several variables. 



Multiple orthogonal polynomials arise naturally in the study of Hermite-Pade approx- 
imation, which is simultaneous rational approximation to a vector of r functions. In 
this respect it is quite natural to study multiple orthogonal polynomials as orthogonal 
vector polynomials. This approach is very useful in trying to extend results for the 
case r = 1 to the case r > 1 by looking for an appropriate formulation using vector 
algebra. Van Iseghem already used this approach to formula a vector QD-algorithm for 

Several algebraic aspects of multiple orthogonal 

A further general- 



multiple orthogonal polynomials [|| 



polynomials follow easily from the vector orthogonality 
ization is to study matrix orthogonality, where the matrix need not be a square matrix 
43|. Orthogonal polynomials and Pade approximants are closely related to certain 



continued fractions (J-fractions and S-fractions). For multiple orthogonal polynomi- 
als there is a similar relation with vector continued fractions and the Jacobi-Perron 
algorithm The seven families which we considered in this paper lead to seven 

families of vector continued fractions, which could be studied in more detail in the 
framework of continued fractions. Finally, one may wonder whether it is possible to 
use hypergeometric functions of matrix argument in the study of multiple orthogonal 
polynomials? 



4.2 Non-symmetric banded operators 

In Section 2 the connection between multiple orthogonal polynomials and banded Hessenberg 
operators of the form 



/ 0-0,0 1 
0-1,1 Oifi 1 
02,2 0-2,1 02,0 



\ 



Or,0 1 

Or+1,0 1 



''n,r—l 



On,l On,0 



\ 



was explained. For ordinary orthogonal polynomials the operator is tridiagonal and can 
always be made symmetric, and often it can be extended in a unique way to a self-adjoint 
operator (e.g, when all the coefficients are bounded). The spectrum of this tridiagonal 
operator corresponds to the support of the orthogonality measure, and the spectral measure 
is precisely the orthogonality measure. Each tridiagonal matrix with ones on the upper 
diagonal and positive coefficients on the lower diagonal, corresponds to a system of orthogonal 
polynomials on the real line (Favard's theorem). Some preliminary work on the spectral 
theory of the higher order operators (r > 1) was done by Kalyagin [^, [^], 0, but 
there are still quite a few open problems here. 
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1. What is the proper extension of Favard's theorem for these higher order banded Hes- 
senberg operators? Not every banded Hessenberg operator corresponds to a system of 
multiple orthogonal polynomials with orthogonality relations on the real line. There 
needs to be additional structure, but so far this additional structure is still unknown. 
There is a weak version of the Favard theorem that gives multiple orthogonality with 
respect to linear functionals ([^, |24|), but a stronger version that gives positive mea- 



sures on the real line is needed. How do we recognize an Angelesco system, an AT 
system, or one of the combinations considered in from the recurrence coefficients 



(from the operator)? The special case where all the diagonals are zero, except for the 
upper diagonal (which contains I's) and the lower diagonal, has been studied in detail 
in 0. They show that when the lower diagonal contains positive coefficients, the oper- 
ator corresponds to multiple orthogonal polynomials on an (r -|- l)-star in the complex 
plane. Using a symmetry transformation, similar to the quadratic transformation that 
transforms Hermite polynomials to Laguerre polynomials, this also gives an AT system 
of multiple orthogonal polynomials on [0, oo). 

2. The asymptotic behavior of the recurrence coefficients of the seven systems described 
above is known. Each of the limiting operators deserves to be investigated in more 
detail. The limiting operator for the Jacobi-Pineiro polynomials is a Toeplitz operator, 
and hence can be investigated in more detail. See, e.g., for this case. Some of the 
other limiting operators are block Toeplitz matrices and can be investigated as well. 
Are there any multiple orthogonal polynomials having such recurrence coefficients? 
The Chebyshev polynomials of the second kind have this property when one deals 
with tridiagonal operators. 

3. The next step would be to work out a perturbation theory, where one allows certain 
perturbations of the limiting matrices. Compact perturbations would be the first step, 
trace class perturbations would allow us to give more detailed results. 

4.3 Applications 

1. Hermite-Pade approximation was introduced by Hermite for his proof of the transcen- 
dence of e. More recently it became clear that Apery's proof of the irrationality of ({3) 
relies on an AT system of multiple orthogonal polynomials with weights Wi{x) = 1, 
u]2{x) = — \og{x) and W3(x) = log^(x) on [0, 1]. These multiple orthogonal polynomials 
are basically limiting cases of Jacobi-Pineiro polynomials where ao = = ai = a2- A 
very interesting problem is to prove irrationality of other remarkable constants, such as 
C(5), Catalan's constant, or Euler's constant. Transcendence proofs will even be better. 



See m ^g] for the connection between multiple orthogonal polynomials, irrationality. 



and transcendence. 

2. In numerical analysis one uses orthogonal polynomials when one constructs Gauss 
quadrature. In a similar way one can use multiple orthogonal polynomials to construct 
optimal quadrature formulas for jointly approximating r integrals of the same function 
/ with respect to r weights Wi, . . . ,Wr. See, e.g., Borges who apparently is not 
aware that he is using multiple orthogonal polynomials. Gautschi [|l^ has summa- 
rized some algorithms for computing recurrence coefficients, quadrature nodes (zeros 
of orthogonal polynomials) and quadrature weights (Christoffel numbers) for ordinary 
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Gauss quadrature. A nice problem is to modify these algorithms so that they compute 
recurrence coefficients, zeros of multiple orthogonal polynomials (eigenvalues of banded 
Hessenberg operators) and quadrature weights for simultaneous Gauss quadrature. 
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